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Abstract 



We introduce a continuous-time quantum walk on an ultrametric space correspond- 
ing to the set of p-adic integers and compute its time-averaged probability distribution. 
It is shown that localization occurs for any location of the ultrametric space for the 
walk. This result presents a striking contrast to the classical random walk case. More- 
I over we clarify a difference between the ultrametric space and other graphs, such as 

cycle graph, line, hypercube and complete graph, for the localization of the quantum 
. case. Our quantum walk may be useful for a quantum search algorithm on a tree-like 

! hierarchical structure. 

O 

o 

Ph! 1 Introduction 



Continuous-time classical random walks on ultrametric spaces have been investigated by 



^ . some groups to describe the relaxation processes in complex systems, such as glasses, clusters 
and proteins [III21E1I11E1EIEIIH]- On the other hand, a continuous-time quantum walk, 
which is a quantum version of the classical random walk, has been widely studied by various 
^ ' researchers for some important graphs, such as cycle graph, line, hypercube and complete 
■ graph, with the hope of finding a new quantum algorithmic technique [HI [IDl [HI [iSl [131 [lH 
[T5l [T6l [T71 [Is] . However any result is not known for a quantum walk on the ultrametric 
space. This is a motivation for our present work. Excellent surveys of quantum walks are 
given in Refs. 

Let X be a set and let p be a metric on X. The pair (X, p) is called a metric space. 
Moreover if p has the strong triangle inequality: 

p{x,z) <max{p{x,y),p{y,z)), 

for any x,y,z & X, then p is said to be an ultrametric. A set endowed with an ultrametric is 
called an ultrametric space. Let p be a prime number and Qp be the field of p-adic numbers. 
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Zp denotes the set of p-adic integers. Then Zp is a subring of Qp and Zp = {x G Qp : \x\p < 1}, 
where | • |p is the p-adic absolute value. Remark that the metric Pp{x,y) = \x — y\p is an 
ultrametric. Every x G Zp can be expanded in the following way: 

X = yo + yip + y2P^ H h ?/„j9" H — , 

where ?/„ G {0, 1, . . . ,p — 1} for any n. See Ref. [21] for more details. It is important that 
Zp can be represented by the bottom infinite regular Cayley tree of degree p, T^'p\ in which 
every branch at each level splits into p other branches. So we consider a continuous-time 
quantum walk on the bottom of T^^^ at level (or depth) M, denoted by T^^ . 

We calculate the probability distribution of the quantum walk on T^^ . From this, we 
obtain the time- averaged probability distribution on and then take a limit as M ^ cxd. 
On the other hand, first we take the limit as M — oo and then we derive the time-averaged 
probability distribution. It is shown that the first limit result coincides with the second one, 
that is, time-averaged operation and M-limit one are commutative. By using the results, 
we clarify a difference between classical and quantum cases. A striking difference is that 
localization occurs at any site for a wide class of quantum walks. Furthermore we compare the 
results of ultrametric space with those of other graphs, such as cychc graph, line, hypercube 
and complete graph in the quantum case. Some applications of p-adic analysis in physics 
and biology were reported [21]. Social network models based on the ultrametric distance 
were presented and studied [221 [23]- The disease spreading on a hierarchical metapopulation 
model was investigated [21]. Our quantum walk on the ultrametric space may be useful for 
designing a new quantum search algorithm on the graph with a hierarchical structure. 

This paper is organized as follows. Section 2 is devoted to the definition of a continuous- 
time quantum walk on t|^\ Our results are shown in Sec. 3. We review corresponding 
results for continuous-time classical random walks in Sec. 4. Section 5 contains results on 
other graphs for the quantum case. Conclusion is given in Sec. 6. 

2 Definition of Quantum Walk 

For M > and x G Zp, the closed ]3-adic ball b'^^[x) of radius p^'^^ with center x is defined 

by 

B%\x) = {ye'Lp:pp{x,y)<p~''}. 

Each ball B^^] (x) of radius can be represented as a finite union of disjoint balls -B^+i {^m) 
of radius 

p-i 

(•^) ~ U ^M+l(^m)' 
m=0 

for suitable Xq, Xi, . . . , Xp_i G Zp, so we have 

pM_i 

Zp= U i?S?(m). 

m=0 
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When p — 3, 



3^-1 3^-1 
m=0 k=Q 

and 



Sf)(0) = 5f (0) U ^(3) U (6), 

fif (l)=5f (l)U5f (4)UBf (7), 
5f)(2) = 5^(2) U £f (5) U 5^(8). 

We define the distance between two balls Bi and B2 as 

S2) = mi{pp{x, y) : X e Bi,y e B2}. 

Wc observe that end points of a regular Cayley tree with p degree at level M may be 
represented set of disconnected balls {B'f^ (0) , B'^^ (1) , . . .,B^^{p^ - 1)} with radius 
p~^ covering Zp. So we can consider balls {5^'*(0), . . . ,B^(p^ — 1)} as p"^ 

points {0, 1, . . . ,p^^ — 1}, denoted by 

As in the case of classical random walk, let (i, j) be amplitude of a jump from the ball 
B^^\i) to the ball B^^\j) separated by p-adic distance p'^^"^). Here we put = e^{k) — 

i^Mihj)- Remark that = ipMij,i)- 

We consider p — 3 and M = 2. In this case, is given by 

ei = (0,3)= (0,6) = 4=^) (3,0)= V^f (3,6) 

= (6, 0) = 4f{6, 3) = V^f (1, 4) = 4'\l, 7) = • • • , 

e2 = 4'^(0,l)=4=')(0,4) = .-- = ^f (0,8) 

= (3, 1) = V'?) (3, 4) = • • • = V'f (3, 8) = • • • , 



Let /„ and J„ denote the n x n identity matrix and the all-one nxn matrix. We define a 
pM y. pM symmetric matrix on T^^ which is treated as the Hamiltonian of the quantum 
system as follows: for any M — 1,2,..., 

n'f^ = tolp + ei{Jp — Ip), 

where €j e R (j = 0, 1, 2, . . .) and R is the set of real numbers. 

The evolutii 
unitary matrix: 



The evolution of continuous-time quantum walk on T^^ is governed by the following 



C/S)(t)=e^*<. 
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The amplitude wave function at time t, |^]^ (t)), is defined by 

In this paper we take |\I/^''(0)) = [1, 0, 0, . . . , 0]"^ as an initial state, where T denotes the 
transposed operator. 

The (n + l)-th coordinate of |\E'j^''(t)) is denoted by t)) which is the amplitude 

wave function at site n at time t for n = 0, 1, . . . — 1. The probability finding the walker 
is at site n at time t on T^^f is given by 

Pi^\n,t) = {¥;^\n,t)\¥ff{n,t)). 



3 Our Results 

Let {rjm : m = 0, 1, . . . , M} be defined by 



Co — ei if m = 0, 

Vm=<, eo + (P - 1) Er=iP'"'efe - P'^em+i if m = 1, 2, . . . , M - 1, 

if m = M. 

Direct computation yields that {77^ : = 0, 1, ... , M} is the set of eigenvalues of 
and eigenvectors for 77^ are given in the following way, see Ref. [T] for p = 2 case. Let 
u = Up = exp{27ii / p) . Put 




The p {p — 1) eigenvectors for rjo are 



Uo{l) = [1, u;^ . . . , \ OpM_pf, 
uo{2) = [l,u;\u;\...,u;'^^-'\OpM_pf, 



uoip - 1) = [1, u;^^-'\uj'^P-'\ 0,M_,] 
uoip) = [Op, 1, t^, . . . , ujP-'^, 0pM_2„] 



T 



Mo(p + 1) = [Op, 1, cu^, . . . , OpM^,pf, 



Uq 



(2(p-l)) = [0„l,>-l),c.2(^-l),...,>-^)^0pM_2p]^ 



Uo{p'''\p - 1)) = [OpM_p, . . . ,c.(^-i)T. 
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The p^^ "^{p ~ 1) eigenvectors for rji are 

mil) 

ui{2) 
ui{p - 1) 

Ui{p) 
Ui{p + 1) 

ui{2{p-l)) 
u,{p''-\p-l)) 

We have — 1) eigenvectors for 77^ (A; = 2, 3, . . . , M — 1) similarly. So the p — 1 

eigenvectors for rjM-i are 

MAf-l(l) = [IpM-l, CUlpM-l, (UJ^lpAZ-l, . . . , """IpM-l]"^, 

Ujv/_i(2) = [IpAf-i , cu^lpM-i , cj^lpAf-i , . . . , a;^'-^ "'"■'ipAf-i]"^, 

mm-i(p - 1) = [V-i,u;^^"^hpAf-i,u;2(p-i)y^_^^ _ _ . ,u;(?'-i)'lpM-i]^. 

Finally the eigenvector for rjM is only [IpM]-^. 
Throughout this paper, we assume that 

< Cm < cm-i < ■ ■ • < £2 < ei. 

The assumption implies that our quantum walk class does not belong to a class of continuous- 
time quantum walks given by the adjacency matrix of the graph on which the walk is defined, 
see Ref. [18]. 

Remark that the diagonal component eo is an irrelevant phase factor in the wave evolution. 
So the probability distribution at time t, {P|^^(n, t) : = 0, 1, . . . — 1}, does not depend 
on eo and we put 

M 

6o = -(p-l)^/-^efe<0, 

k=l 

as in the classical case. We should note that eo < < ejv/ < ^m-i < ■ ■ ■ < e2 < ei is 
equivalent to r^o < < " " " < Vm-i < Vm = 0. 

Let V^^'^ = {p'^-i,/-! + 1, . . . , / - 1} (A; = 1, 2, . . . , M) and V^^^^ = {0}. Remark that 
Ujl^ Vfc(p) = {0, 1, . . . - 1} and = {p- for A: = 1, 2, . . . , M, where \A\ is the 

number of elements in a set A. By using the eigenvectors, we obtain the amplitude of the 
quantum walk on T^"*: 



[lp,ujlp,uj Ip, . . . ,uj^ Ip, Opj\/_.p2] , 

[Ip, Uj'^lp, Uj'^lp, . . . , Uj'^^^'^hp, OpM_p2]^ 



[lp,uj^P-'hp,co'^P-'hp, . . . ,u;(P-i)'lp,0pM_p2]^, 
[0p2, Ip, iolp, iLp'lp, . . . , u;^~"^lp, OpA/_2p2]^, 

[0p2, Ip, iO^lp, LO^lp, . . . , (^^^^"""^^Ip, 0pM_2p2]^, 

[0p2, Ip, uo^^-^hp, uo'^'^^'^hp, . . . , uo^P-^^hp, OpM_p2]^, 



[OpM_p2, lp,a;(^'-i)lp,a;2(p-i)i^, . . . ,o;(P-i)'lp]^ 
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Lemma 3.1 



M-l 



■M 



m=0 



ip) 



M 



m=k 



^-M^_^rtr,M-r + 1) 



tfneVif^ (fc = l,2,...,M-l), 
) 



The definition of P^\n,t) implies 
Proposition 3.2 



Af-l 



m=0 



M 



'M-l 



p 



-(m+l) 



sin(t?7„ 



(p) 



. m=0 



M-l 



-p-'' cos(tr/fc-i) + {p-l)J2 cos(tr7™) +p 



m=k 



AI-1 



^/r^GFi^) (A; = 1, 2, . . . , M - 1), 



m=k 



2p-^^^(l-cos(tr^,,_J) 
>(p)/ 



If n = 0, then P^''(0,t) is the return probability of the walk on t|,j\ For p = 3 and 
M = 2 case, we obtain 

' {41 + 24 cos(t(r/o - r?i)) + 12 cos(tr7o) + 4 cos(tr7i)} /3^ 
if n = 0, 

{14 - 12cos(t(r7o - Vi)) -Qcos{tr]o) + 4cos(tr/i)} /3^ 

if n = 1,2, 
2{l-cos(tr/i)}/3^ 

if n = 3,4, 



Pf (n,t) = <( 



^, J., . . . , . 



In general, P^\n,t) does not converge as t ^ oo for any fixed n. So we introduce a 
time-averaged distribution of P^j{n,t) as follows: 

1 ''^ 



if the right-hand side of Eq. (13.11) exists. Then Proposition 13.21 gives 



>(p)/ 



(3.1) 
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Theorem 3.3 



p-1 2 

1 1 

+ 



,2M 



ifneVlf^ (A; = 1,2,...,M-1) 



It is interesting to note that P^\n) does not depend on {ek : /c = 0, 1, . . . , M}. In the case 
of p = 3 and M = 2, we have 



/^f (0) = |, /-f (1) = (2) = ^, Pf (3) 
The following result is immediate from Theorem 13.31 
Corollary 3.4 
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2 

34- 



hm P^\n) 

M^oo 



r{p) 



ifneV^"' (A: = l,2,...). 



{p + l)p2fc-l 

Here we consider the mean distance from at time t defined by 

M 



k=l 



(p) 



Then its time-averaged mean distance is given by 



M 



n] 



k=l 



(p) 



From Theorem 13.31 we get 

_ 2(j9-l)(M-l) , 2[{(p-l)(p+l)^ + l}p^ 



2M-2 



d 



1] 



M 



(p + lyp 



,3M-1 



This gives 



M^oo M p+1 



Next we take a limit as M ^ 00 first. We define Pi?'' (n, t) by Pi?'' (n, t) = limA/^00 Pm'' {n, t) , 
if the right-hand side of the equation exists. By Proposition 13. 2^ we obtain 
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Proposition 3.5 



P^\n,t) = { 



.m=0 



P 



-(m+1) 



cos(trir, 



(m+1) 



sin(t?7n 



, m=0 



cos(t?7fc-i) + (p - 1) ^ cos(i(;r7„ 



+ <^ -p-' sin(%_i) + (p - 1) ^ sin(tr7„ 

^. m=k 

ifneVf^ (A: = l,2,...). 
In a similar way, a time-averaged distribution P,^'' (n, t) is given by 

Pi^)H= liml fp!^\n,s)ds, 

t^oo t Jq 



(3.2) 



if the right-hand side of Eq. (13.21) exists. Combining Theorem 13 . 31 with Proposition 13 .51 yields 
Corollary 3.6 For any n = 0,1,2, .. ., we have 

Furthermore, 

\imP(^\n) = 6o{n), 



lim Pi?(n) > 0. 

M^oo 



p—*oo 



where Sm{n) = 1, if n = m, = 0, if n 7^ m. 

In general, we say that localization occurs at site n if the time-averaged probability at 
the site is positive. Therefore the localization occurs at any site for both any finite M and 
M —>■ 00 limit cases. 



4 Classical Case 

In this section we review three classical examples and clarify a difference between classical 
and quantum walks. Let Pjf^ (n, t) be the probability that a classical random walker starting 
from is located at site n at time t on Zp. 

In the case of a linear landscape, the transition rate on T^^ has the form 
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for wq > and a > 0. Then Avetisov et al. [¥] showed a power decay law in M — ^ oo hmit: 

^i^nO,t)~^ (t-oo), 
where f{t) ~ g{t) {t oo) means there exist positive constants Ci and C2 such that 

Ci < hminf 44 < hmsup 44 < C's- 

i^oo g[t) t^oo g[t) 

For a logarithmic landscape case, the transition rate on T^"* is 

eu = Wop-^'-''^ ^ 



for i(7o > and a > 1. The following stretched exponential decay law (the Kohlrausch- 
Williams- Watts law) was proved by Avetisov et al. |5j in M — >■ cx) limit: 

log(PiP)(0,t))~-ti/" (t-.oo). 

In the case of an exponential landscape, the transition rate on T^^^ has the form 

= «;op-('=-A^)exp(-ap('=-^'^)), 

for > and a > 0. Then Avetisov et al. [5j obtained a logarithmic decay law taking a 
limit as M — > 00: 

i^i^nO,t)~r^ (t-00). 
logt 

The above three facts imply that the localization does not occur at position 0. In contrast 
to the classical case, the localization occurs at any position for our quantum case. 



5 Quantum Case for Other Graphs 

We consider the time-averaged probability distribution for continuous-time quantum walk 
starting from a site on other graphs, such as cycle graph, line, hypercube and complete 
graph. Then the Hamiltonian of the walk is given by the adjacency matrix of the graph. 
In the case of a cycle graph with N sites, we have obtained the following result 



N m 



for any n = 0, 1, . . . , — 1, where 



if AT = odd, 6n7^0 (mod27r), 
Rj^(n) = { -I a N = even, ^2n ^ (mod27r), 

if6n = (mod27r). 
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Here C,j = 2TTj/N, N = [{N — l)/2], and [x] is the smallest integer greater than x. When 
N = odd (i.e., N = {N-l)/2), 



N-l 



1 1 1\ 



when N = even (i.e., = (A^ - 2)/2), 



(N-2)/2 (7V-2)/2 



Z : ^ . . 'T ^ 



1 12 12 1 A^-2 12 ^ ^ \ 

Then we have 



lim Pjv(^) = 0, 

Af— >oo 

for any n. 

In Z case, the probability distribution P{n, t) of the walk starting from the origin at 
location n and time t is given by 

P(n,t) = J2(t), 

where J„(t) is the Bessel function of the first kind of order n, see Ref. [T3], for example. The 
asymptotic behavior of J„(t) at infinity is as follows: 

fY / An^ - I \ 

Jn{t) = d—icos{t-e{n))-sm{t-e{n))—^ + 0{t-^)\ , t ^ oo, 

where 6{n) = {2n + l)7r/4, see page 195 in Ref. |25j. From this fact and | Jn(^)| < 1 for any 
t and n (see page 31 in Ref. [2S]), we obtain 

1 f 



1 /-^ 

P{n) = lim - / J^(s) rfs = 0, 

t^OO t Jq 



for any n G Z. 

Next we consider the quantum walk starting from a site, denoted by 0, on a hypercube 
with 2^ sites, W^- Then the probability finding the walker at site n at time t is 

P^(n, t) = cos(t/A^)2(^-^) sin(t/A^)2'= if n G 14 (0 < A; < A^), 

where = {x G PVat : | |x| | = k} and | |a;| | is the path length from to a; in Wn. A derivation 
is shown in Ref. Il8i . From the result we see that 



- , , 1 fN\ f2k\ f2(N - k) 



2'^^\kJ \k J \ N -k 

if n G Vfc (0 < /c < A^). For example, when A^ = 4, 

35/128 if n G I/q U V4, 
P^{n) = { 5/128 if n G ^1 U V^, 

3/128 if n G V2. 
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In general A^, we have 



lim P7v(n) = 0, 

Af— >oo 



for any n = 0, 1, . . .. Moreover if we let 



neVfe 



for k = 0,1, . . . , N, then we get 



^^\k J\ N-k 

since \Vk\ = N\/{N — k)\k\. It is interesting to note that this probability corresponds to a 
well-known arcsin law for the classical random walk. 

Therefore, in three cases, Cat (cycle graph) and Wn (A^-cube) as A^ — >• oo and Z, the 
localization does not occur at any location. 

In the case of the complete graph with A^ sites, K^, the Hamiltonian if^v is defined by 
if TV = In — NJ]\r. Then eigenvalues of if at are rjo = 0, rji = ri2 = ■ ■ ■ = rj^_i = —N. The 
eigenvectors are given by the Vandermonde matrix, see Ref. [H]. Then direct computation 
yields 

(AT - 1)2 + 1 + 2{N - 1) cos{Nt) 

if n = U, 



PN{n,t) = { 2(l-cos(Art))^' 



iV2 



a n = l,2,...,N -1. 



Therefore 



(AT 



if n = 0, 



if n = 1,2, ...,A^- 1. 



So we have 



lim PNin) = 5o{n), 

N^oo 



for any n = 0, 1, . . . . Therefore the localization occurs at only n = site. This corresponds 
to our case for p — cxd. 



6 Conclusion 

We have derived the expression of the probability distribution of a continuous-time quantum 
walk on t|^'' corresponding to Zp in the M — >■ cxd limit. As a result, we obtained 

P^\n)= lim Pi?(n)>0, (n = 0, 1, 2, . . .), 
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for a class of satisfying 

to < < tM < CM-l < ■ • • < £2 < ei. 

Therefore the locahzation occurs at any location. For Cn (cycle graph) and Wn (A^-cube) 
as — > oo and Z cases, the localization does not happen at any location. For (complete 
graph) as ^ oo case, the localization occurs at only site. In three typical classical cases, 
the localization does not occur even at site. We hope that this property of our quantum 
walk can be useful in a search problem on a tree-like hierarchical structure. 
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